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We study the problem of joint estimation of real squeezing and amplitude of the radiation field, 
deriving the measurement that maximizes the probability density of detecting the true value of the 
unknown parameters. More generally, we provide a solution for the problem of estimating the un- 
known unitary action of a nonunimodular group in the maximum likelihood approach. Remarkably, 
in this case the optimal measurements do not coincide with the so called square-root measurements. 
In the case of squeezing and displacement we analyze in detail the sensitivity of estimation for coher- 
ent states and displaced squeezed states, deriving the asymptotic relation between the uncertainties 
in the joint estimation and the corresponding uncertainties in the optimal separate measurements of 
squeezing and displacement. A two-mode setup is also analyzed, showing how entanglement between 
£^ , optical modes can be used to approximate perfect estimation. 
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> ' I. INTRODUCTION 

m ' 

, Squeezing and displacement are the basic operations of continuous variables quantum information 1] , and are easily 
performed, the former by parametric amplifiers, the latter by lasers and linear optics. Squeezing in the two-mode 
setup is, for example, the tool to generate entanglement in the Braunstein-Kimble teleportation scheme Q. The 
combined use of real squeezing and displacement allows one to encode efficiently classical information in quantum 
channels using homodyne detection at the receiver 3 J . In a quantum communication scenario where a coherent signal 
is sent through a non-linear medium and undergoes an amplification process, the joint estimation of displacement and 
O |' squeezing provides a twofold information: the amplitude modulation of the state and a property of the communication 
channel itself. This may be useful in a communication scheme designed to be robust to photon loss. 

In this paper we consider the problem of jointly estimating real amplitude and real squeezing of the radiation 
field. In a tomographic setup, where a large number of equally prepared copies is available, the maximum likelihood 
^ method turns out to be very efficient in estimating the parameters that characterize the state of the quantum system 
0, . In this approach, first one fixes a set of single-copy measurements (typically homodyne measurements at some 
random phase) and then looks for the estimate that maximizes the probability (density) of producing the observed 
data. However, the tomographic approach is not suitable to the case when only a small number of copies is available, 
and one needs to use the limited resources at disposal more efficiently. It becomes then important to optimize 
not only the posterior processing of the experimental data, but also the choice of the measurement that is used to 
extract these data from the system. The most natural framework to deal with this situation is quantum estimation 
theory ||J Q , where the concept of positive operator valued measure (POVM) provides a tool to describe at the same 
time both measurement and data processing. The maximum likelihood approach in quantum estimation theory Q 
then corresponds to seek the measurement that maximizes the probability (density) that the estimated value of the 
parameters coincides with the true value. 

Joint estimation of squeezing and displacement is equivalent to infer an unknown transformation of a group — in 
the present case the affine group. This is an example of a frequent situation in quantum estimation, especially in 
communication problems, where a set of signal states is generated from a fixed input state by the action of a group. 
Consider, for example, the case of phase estimation for high-sensitivity interferometry and optimal clocks P. llOllTl| . 
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the estimation of rotation for the optimal alignment of reference frames [12|, and the estimation of displacement 
of the radiation field for the detection of a coherent signal in Gaussian noise 0]. In these cases, the symmetry 
of the problem provides a physical insight that allows one to simplify the search for efficient estimation strategies, 
and the concept of covariant measurement becomes crucial for optimization. The use of maximum likelihood 
approach in the covariant setting has been shown to be particularly successful in obtaining explicitly the optimal 
measurements and in understanding the fundamental mechanism that leads to the ultimate sensitivity of quantum 
measurements 0,0]. Group theoretical tools, such as equivalent representations [l4| and multiplicity spaces [l5| . 
far from being abstract technicalities, are the main ingredients to achieve the ultimate quantum limits for sensitivity 
|l2l Il6| . Moreover, in a large number of situations, the maximization of the likelihood provides measurements that 
are optimal also according to a wide class of different figures of merit [16| . 

From a group theoretical point of view, the action of real squeezing and displacement on the wavefunctions provides 
a unitary representation of the affine group "ax + 6" of dilations and translations on the real line. The structure of 
the affine group underlies the theory of wavelets ^tJj an d has been recently used in the characterization of SU(1, 1) 
coherent states 0, , and in the study of oscillators in a Morse potential (2j| • 

The affine group is particularly interesting, since it is the paradigmatic example of a nonunimodular group, namely 
a group where the left-invariant Haar measure is different from the right-invariant one. This leads to orthogonality 
relations that — differently from the usual Schur lemmas — involve a positive unbounded operator, firstly introduced by 
Dufio, Moore, and Carey (DMC) in Ref. 21] . As we will see in this paper, the nonunimodularity of the group has some 
amazing consequences in the estimation problem. For example, the so-called square-root measurements 22], that are 
commonly considered in quantum communication and cryptography, do not coincide with the maximum-likelihood 
measurements, the latter providing a higher probability density of correct estimation. Another bizarre feature is that 
for maximum likelihood measurements, the most likely value in the probability distribution can be different from the 
true one. While for unimodular groups this feature never happens, for nonunimodular groups it is unavoidable, and 
a suitable choice of the input states is needed to reduce the discrepancy between the true value and the most likely 
one. 

The paper is organized as follows. In order to set the optimal joint estimation of real squeezing and displacement in 
the general estimation method, first we derive in Sec. [H]the optimal measurement for the problem of estimating the 
unknown unitary action of a nonunimodular group. Then, as a special example, the general results will be used to 
optimize the joint estimation of squeezing and displacement in Sec. 11111 The efficiency of coherent states and displaced 
squeezed states is analyzed in detail, and the asymptotic relation between the uncertainties in the joint estimation 
and the uncertainties in the optimal separate measuremnts of squeezing and displacement is derived. The conclusions 
are summarized in Sec. IV. The explicit derivation of group average over the affine group is given in the Appendix. 

II. MAXIMUM LIKELIHOOD ESTIMATION FOR A NONUNIMODULAR GROUP 

A. Background and generalities 

Suppose that a fixed input state p, corresponding to a density operator in the Hilbert space 7i, is transformed by 
the unitary representation {U g \ g € G} of the group G, so that it generates the family of signal states 

O = {U gP Ul I g e G} . (1) 

The typical quantum estimation problem is then to find the measurement that gives the best estimate for the unknown 
transformation g according to some optimality criterion. Usually the criterion is given by a cost function c(g,g), 
which quantifies the cost of estimating g when the true value is g, and enjoys the invariance property c(hg, hg) — 
c (<?jS) V/i, <7, g € G. For example, the maximum likelihood criterion corresponds to the delta cost-function c(g,g) = 
—S(g~ 1 g) (loosely speaking, there is an infinite gain if the estimate coincides with the true value, and no gain 
otherwise). Once a cost function is fixed, one can choose two possible approaches to optimization, namely the Bayes 
approach and the minimax. In the Bayes approach, one assumes some prior distribution of the unknown parameters, 
and then minimizes the average over the true values of the expected cost c(g) = J dg c(g, g)p{g\g), where p(g\g) is the 
conditional probability density of estimating g when the true value is g. In the minimax approach, one looks instead 
for the measurement that minimizes the supremum of the expected cost over all possible true values of the unknown 
parameters. 

An important class of estimation strategies is given by the covariant measurements, that are described by POVMs 
of the form 



M(g) = UgZUl , 



(2) 
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where £ > is an operator satisfying the normalization condition 

/ d L gU g €Ut = t, (3) 
Jg 

d-L g denoting the left-invariant Haar measure on the group, namely dL(hg) = d^ g V7i, jeG. Due to the symmetry 
of the set of states Q , covariant measurements play a fundamental role in the search of the optimal estimation. For 
compact groups, the following proposition holds: 

Proposition 1 (Holevo 0)- for compact groups, the search j or the optimal measurement in the minimax approach 
can be restricted without loss of generality to the class of covariant measurements. 

Moreover, for compact groups the optimality of covariant measurements holds also in the Bayes approach 0, if 
the prior distribution is chosen to be the normalized Haar measure on the group, i.e. the measure dg such that 
dg = d{hg) = d{gh) Wg,heG and J G dg = 1. 

For non-compact groups, such as the affine group involved in the joint estimation of squeezing and displacement, 
the situation is more involved. Of course, in the Bayes approach it is no longer possible to choose the uniform Haar 
measure as prior distribution, since it is not normalizable. However, in the minimax the optimality of covariant 
measurements still holds, even though in this case the proof becomes rather technical j^. In this paper we will adopt 
the minimax approach to deal with noncompact groups, and this will allow us to restrict the optimization to covariant 
POVMs. 



B. Nonunimodular groups 

As we will see in Section III, the action of real squeezing and displacement on one mode of the radiation field 
yields a representation of the affine group u ax + 6" of dilations and translations on the real line. This group is clearly 
non-compact, and, moreover, it is nonunimodular, namely the left-invariant Haar measure dr,g (with di,(hg) — 
d^g Vh,g G G) does not coincide with the right-invariant one dug (with du(gh) = dug Vh,g G G). Therefore, to 
face the estimation problem with the affine group, we need some results about representation theory and orthogonality 
relations for nonunimodular groups (for an introduction to these topics, see for example |24|V 

Let {Ug | g G G} be a unitary representation of a locally compact group G. In the following, we will make two 
assumptions on the representation {U g } that are tailored on the concrete problem of estimating real squeezing and 
displacement. 

First assumption: discrete Clebsch-Gordan series. We require the representation {U g } to be a direct sum of 
irreducible representations (irreps), namely that its Clebsch-Gordan series is discrete. In this case, there is a decom- 
position of the Hilbert space H as 

« = Hi* ® M„ , (4) 
pes 

such that 

Vg = u 9 ® > (5) 
pes 

where {Ug 1 } is an irreducible representation acting on the Hilbert space H M and 1^ is the identity in the space M.^. 
The Hilbert spaces 7i M and M.^ are called representation and multiplicity spaces, respectively, and the index /i labels 
the inequivalent irreps that appear in the Clebsch-Gordan series S of the representation {U g }. 

Second assumption: square- summable irreps. We require each irreducible representation {U^} in Eq. JSJ to be 
square-summable. This means that there is at least one non-zero vector |^) G TL^ such that 

/ d i9 |<Vv|^|VV>| 2 <°°. (6) 
Jg 

Vectors such that the above integral converges are called admissible. It is possible to show [24| that, if a representation 
is square-summable, then the set of admissible vectors is dense in the Hilbert space. 
Let us consider now the group average of an operator A, defined as 



(7) 



4 



In general, the average may not converge for any operator (for example, it diverges for A = 1). In analogy with 
admissible vectors, we say that A is an admissible operator if the group average in Eq. converges in the weak 
operator sense. In such a case, one can prove that (A)g is given by 

(A) G = In, ® Tr Hlt [D„ ® l Mtt A] , (8) 

where D M is a positive self-adjoint operator actin g on the representation space TL^. The operator has been firstly 
introduced by Duflo, Moore, and Carey (DMC) [21|. and is the characteristic feature of nonunimodular groups. In 
fact, if the group G is unimodular — i.e. if the left- and right-invariant measures coincide — than the DMC operator 
is simply a multiple of the identity, and the formula © for the group average is equivalent to the ordinary Schur 
lemmas. Contrarily, if the group G is nonunimodular, the DMC operator is a positive unbounded operator, and 
its presence modifies the orthogonality relations dramatically, with remarkable consequences in the estimation of an 
unknown group transformation. 

The admissibility of vectors and operators has a simple characterization in terms of the DMC operator. As regards 
vectors, the set of admissible vectors for the irrep {U^} can be characterized as the domain of y/D^. For unimodular 
groups, since is proportional to the identity, the set of admissible vectors is the whole representation space Ti^, 
while for nonunimodular groups the admissible vectors form a dense subset of TL^. As regards operators, an operator 
A is admissible if and only if the partial traces Tr-^ [D^ £g) 1m A] are not diverging. For example, if \ipu) G Ti.^ is an 
admissible vector and is any operator acting on A^ M , then the operator A^ = IVvXVvl ® ls admissible, and its 
group average is (A^a = (tp^D^tp^) 1 M ® M . 



C. Maximum likelihood measurements 



In order to find the best estimate for the signal states (JIJ according to the maximum likelihood criterion in the 
minimax approach, we consider now a covariant measurement, described by a POVM M{g) as in Eq. @. The 
normalization condition @ can be rewritten as 

f d L g M(g) = (Z) G = l . (9) 
Jg 

Using Eq. |(SJ| for evaluating the group average (£)gj this condition becomes 

Tr Hit [D^l Mit ^]=t Mit . (10) 

According to the maximum likelihood approach, we need to find the covariant measurement that maximizes the 
probability density that the estimated transformation coincides with the true one, i.e. C = p(g\g) = Tr[M(g)U g pU^] = 
Tr[£p]. For a pure input state p — \ip)(ijj\, the maximization of the likelihood over all possible operators £ > satisfying 
the constraints (|10f) follows in a simple way by a repeated use of Schwartz inequality. In fact, the input state \ip) can 
be written in the decomposition Q as 

IV>>=® c p |* p », (11) 
where each |^)) G Ti^ ® M.^ is a bipartite state. From Schwartz inequality, we have 



fi,V \ fJ, 
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/':,:>. IvvWtl*, ! (>>,K * ! (12) 

1/2 

At this point, we assume that each bipartite state l^^)) is in the domain of the operator 1 ® 1m, ■ This assumption 
is not restrictive, since the domain of a self-adjoint operator is dense in the Hilbert space. In this way, it is possible 
to write |^)) = D^D^ 1 ^ 2 ® IxJ'IV)) and to exploit the Schmidt decomposition of the (non-normalized) vector 

— 1/2 

Dp (g) 1m, I^m))- I n other words, we can write 

= E Vd;® iA4jc)ic) , (is) 

m— 1 
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where r p is the Schmidt rank, A^ n > are Schmidt coefficients such that J2m=i = ((^Vl-^/x 1 ® ^-M J an d 
IV'm)) IVm) are the elements of two orthonormal bases for Ti^ and .M^, respectively. The form ((13(1 is very convenient 
for optimization, in fact we can use again Schwartz inequality and obtain 



< ( E v A ™ $&K&IV^® i^J^)lfc)J (14) 

Finally, we have 

(^kci^®im^vs;®imJC)ic)<(ci twua.® 1 ^] ic) = i, (15) 

the last equality following from the normalization constraint 111) II . Therefore, the previous chain of inequalities proves 
the upper bound 

C< (5>„| £ v^J = (16) 
y^GS m=l / 

that holds for any covariant POVM. On the other hand, it is immediate to check that the bound is achieved by the 
covariant POVM given by £ = |?7)(?7| with 

\V) =0e»^ J2 ® |<)|C> , (17) 

^GS m=l 

arg(z) denoting the argument of a complex number, i.e. z — |z|e larg( ^. The normalization of such a POVM follows 
from Eq. I|l()[). and one has 

J d L 5 C/ 9 £C/t = 8 £ , (18) 

m— 1 

namely, the POVM is complete in the subspace spanned by the orbit of \ip), and can be trivially completed to the 
whole Hilbert space without affecting the probability distribution. Notice that, if the group G is unimodular, namely 
Dp = 1« /dfj, for some positive constant d^, then we correctly retrieve the results of Ref. 15] about maximum- 
likelihood measurements. 

The case of {U g } being a direct sum of inequivalent irreps. The expression for the optimal covariant POVM can be 
further simplified in the case when all the multiplicity spaces M.^ are one-dimensional, i.e. when the representation 
{U g } is a direct sum of inequivalent irreps. In this case, we can decompose the input state as \ip) = c^ip^) (as 
in Eq. 1)11(1 . but without the need of introducing bipartite states), and now the decomposition ((13(1 becomes trivial, 
namely 



IVV) = v -V V D M > (19) 

where A M = (Vvl-D^lVv) and 

= ■ (20) 
n^ i/2 i^)n 

Therefore, Eq. ((17(1 for the optimal POVM becomes 

\„) = e i««(e,) D "^ ] = , (21) 

and the corresponding optimal likelihood is given by 



c°^ = \m\ 2 = [£vw^> I • ( 22 ) 
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D. Remarks 



Remark 1: Square-root measurements. A possible strategy to estimate an unknown quantum state, randomly drawn 
from a given family, is given by the so-called square-root measurements (SRM), firstly introduced by Hausladen and 
Wootters (2^1 - In the case of pure states with a group symmetry as in Eq. Q, there is an important connection 
between SRM and maximum-likelihood measurements. For example, if the group G is a discrete group of phase 
shifts, it has been proved in Refs. [25l I26L . that SRM minimize the probability of error in estimating the unknown 
state, and, equivalently, they maximize the likelihood, i.e. the probability of correct estimation. More generally, the 
optimality of the SRM in the maximum likelihood approach has been proved in Ref. for a large class of groups, 
including all finite groups, all compact groups, and unimodular noncompact groups, such as the Weyl-Heisenberg 
group of displacements. However, as we will see in the following, the case of nonunimodular groups represents an 
exception to the fact that SRM are optimal for the maximum likelihood criterion in the presence of a physical 
symmetry. In fact, the SRM for the estimation of a group transformation acting on a fixed state p is given by the 
POVM M sq (g) = F- 1 ' 2 U g pW g F- 1 / 2 , where 



F = f A L gU gP Ul (23) 

JG 



(the POVM M sq (g) is obviously normalized with respect to the left-invariant measure d^ g). The comparison with the 
maximum-likelihood measurements of the previous section is particularly simple in the case of group representations 
{U g } that are direct sum of inequivalent irreps. In fact, for a pure state p = \^)(tjj\ with 1^) = © M c^ip^), the integral 
(|2*3l is easily calculated by using Eq. (|SJ), namely F — 0^ \Cfj,\ 2 (tp^D^^) ln„- Notice here that the square-root 

1 /2 

measurement can be defined only if is in the domain of Z? M . Therefore, the square-root measurement is given 
by the covariant POVM M sq (g) = U g \r] sq )(r] sq \U^, where 

This covariant POVM is different from the optimal one given in (|21|l . and does not achieve the optimal value l|22|l for 
the likelihood. When |^) is in the domain of both D]J 2 and 1 , we can compare the values of the likelihood as 
follows. One has 



where we used the inequality {tp^D^ 1 ^^) (ip ^\D ^\ip M ) > 1, from Schwartz inequality applied to the vectors Dj/ 2 ^^} 
and D^ 1/2 \^^). 

Remark 2: The true value and the most likely one. In the maximum likelihood approach, one optimizes the choice 
of the POVM in order to maximize the probability density that the estimated value of the parameters coincides with 
the true one. Intuitively, one could expect that the probability distribution p(g\g) = Ti[M(g)U g pU^] for the optimal 
POVM achieves its maximum at the value g — g. Again, this is true for unimodular groups, but fails to hold for 
nonunimodular groups. 

Proposition 2. Let the group G be unimodular. If the covariant POVM M(g) maximizes the likelihood for a given 
input state p, then the probability distribution p(g\g) of the estimate g on the state U g pU^ achieves its maximum for 
9 = 9- 

Proof. Suppose that the most likely value does not coincide with the true one. Then we can rigidly shift the whole 
probability distribution with a post-processing operation that brings the most likely value to the true one. In fact, if 
the maximum of p(g\g) occurs at g — gh, we can always replace M(g) with a new covariant POVM M'{g) = U g £' ut , 
where 

£' = U h £Ul . (26) 

The normalization of the new POVM follows from the fact that for unimodular groups the DMC operators are trivially 
proportional to the identity, and therefore the operator £' satisfy the normalization constraints (|10|l as well. Moreover, 
the probability distribution p'(g\g) associated with M'(g) enjoys the property p'(g\g) = p(gh\g), whence it achieves 
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the maximum in g = g. In this way, the likelihood of M'(g) would be higher than the likelihood of the POVM M(g). 
But this cannot happen since M(g) is the optimal maximum-likelihood POVM. Therefore p{g\g) must be maximum 
in g = g. ■ 

For nonunimodular groups the previous argument does not apply, since the POVM given by (|26|) is no longer 
normalized. In fact, the operator £' does not satisfy the normalization constraints (|10f> . since the DMC operators do 
not commute with the unitaries Uh- In other words, we are not allowed to bring the most likely value to coincide 
with the true one by rigidly shifting the whole probability distribution. As we will see in the explicit example of the 
estimation of real squeezing and displacement, this situation can indeed happen. In order to reduce the discrepancy 
between the true value and the most likely one, a suitable choice of the input states is needed. For example, in the 
simple case of {U g } being a direct sum of inequivalent irreps, if the projection of the input state onto the irreducible 
subspaces are eigenvectors of the DMC operators, then the most likely value coincides with the true one. In fact, for 
any input state \ip) = @ M c m |Vv)j using Schwarz inequality, we have 



p opt (g\g) 



< 



(27) 



and if each | ip^ 
one. 



is eigenvector of D M , then the last expression is equal to p(g\g), then the true value is the most likely 



III. OPTIMAL ESTIMATION OF REAL SQUEEZING AND DISPLACEMENT 

A. Translation and dilation 



In the following we will apply the general framework of Section[n]to the case of joint estimation of real squeezing and 
displacement of a single- mode radiation field with bosonic operators a and a' with [a, at] = 1. Given the wavefunction 
of a pure state \ip) in the A-representation ip( x ) — (^l^); where \x) denotes the Dirac- normalized eigenstate of the 
quadrature operator X = (a + a')/2, the affine transformation on the real line given by x — » e r x + x' is represented 
by the unitary transformation 



i/)(x) er r/2 ijj(e- r (x - x')) , x',r € R . 
This transformation corresponds to the action of the unitary operator U x i ir = D(x')S(r) on the ket where 

D{x) = exp \x(a' — a)] , 



(28) 



S(r) 



exp 



2> t2 



a 2 ) 



(29) 



represent the displacement and the squeezing operator with real argument, respectively. In other words, the operators 



{U x , r = D{x)S{r) I x ,reR} , (30) 

provide a unitary representation of the affine group in the Hilbert space of wavef unctions. The affine group is 
nonunimodular, and in the above parametrization the left- and right-invariant measures are given by d^ g = e~ r d r d x, 
and d^ g = d r d x, respectively. 



B. The maximum likelihood POVM 



In order to exploit the general results of Sec. [HI we need to know the Clebsch-Gordan decomposition of the 
representation {U x , r }, the irreducible subspaces, and the DMC operators. All these informations are given in the 
following, while the proof is presented in the Appendix. 

The Clebsch-Gordan series of the representation {U x ,r} consists on two irreps, that we indicate with the symbols 
+ and — . Accordingly, the Hilbert space splits into two irreducible subspaces, i.e. 



H = H+®H~ 



(31) 
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Comparing this decomposition with the general case J2J, we see that the subspaces Tt+ and Tt- are the representation 
spaces, while the multiplicity spaces M + and M._ are trivially one-dimensional. The representation spaces 7i + and 

TL- can be easily characterized in terms of the quadrature Y = . In fact, writing the wavefunctions in the 
^-representation as ip(y) = (y\ip), where \y) are the Dirac-normalized eigenvectors of Y , we have 7i + = {\ip) | ip(y) 
Vy < 0} and H- = {\tp} \ tp(y) = Vy > 0}. Therefore, the projection operators onto H+ and H- can be written 
respectively as 1+ = 9{Y) and 1_ = 9{— Y), where 9(x) is the customary step-function \9{x) = 1 for x > 0, 9(x) = 
for x < 0]. Moreover, the DMC operators are 

D±=tt . (32) 

With these tools we are now able to provide the optimal covariant measurement for the joint estimation of real 
squeezing and displacement on a given state of the radiation field. Let us denote by \ip) the input state that undergoes 
to unknown squeezing and displacement transformations. Decomposing the input state on the subspaces 7i+ and 7Y_ 
as = c+\ip+)+C-\ip-), we can exploit Eq. it2"T|) and write explicitly the optimal POVM as M(x, r) = U x ^ r \r]}(r]\U^ r 
where 

= \Y\0(YM) \Y\9(-YM) 

^ mmxW) (m\o(-YM ' 

The optimal likelihood is then 

= i (V m \Y\9{Y) W) + vWM¥) 2 . (34) 

according to the general expression of Eq. I|22l) . As already mentioned, the expression of the likelihood provides 
some insight about the states that are most sensitive in detecting an unknown combination of real squeezing and 
displacement. Essentially, one can improve the likelihood by increasing the expectation value of \Y\, the modulus 
of the quadrature Y. In addiction, the use of wavefunctions that in the ^-representation are non-zero both in the 
positive half- line and in the negative half- line allows to exploit the interference of the components \ip+) and \ip-) to 
enhance the value of the likelihood. 



C. Real squeezing and displacement on a coherent state. 

Using Eq. (|33|l for the optimal POVM, we can obtain the probability distribution of the estimated squeezing 
and displacement parameter for a given input state. In particular, for a coherent input state \a) the sensitivity of 
the measurement can be significantly improved by increasing the imaginary part of a, this corresponding to taking 
coherent states with a high expectation value of |V|. The probability distribution for the joint estimation of squeezing 
and displacement on the vacuum state and on a coherent state with a = lOi has been reported in Figs. 1 and 
2, respectively. Comparing the two figures, a remarkable improvement in the precision of the measurement can be 
observed as an enhancement of the likelihood, along with a narrowing of the probability distribution. Moreover, we 
can observe that for the vacuum state the maximum of the probability density is not achieved by the true value (which 
is given by to x = r = 0). The discrepancy between the most likely value and the true one, due to the fact that the 
affine group is nonunimodular, essentially disappears by increasing the expectation value of \Y\. As we can see from 
Fig. 2, for a = lOi the probability distribution is approximately a Gaussian centered around the true value x = r = 0. 

Now it is interesting to focus on the asymptotic behavior of the probability distribution for a coherent state \ia) 
with ael going to infinity. In the asymptotic regime, the probability distribution 

p a (x, r) d x dre~ r = \(r/\ U x ^ r \ia}\ 2 dx dre~ r , (35) 

given by optimal vector \rf) in Eq. I|33|l . can be further simplified. In fact, the wavefunction of the coherent state can 
be written as 

(V\ia) = 0) 1 e- ( ^ Q)2 , (36) 

and for a 3> 1 it is essentially confined in the positive part of the y-axis. Hence, in the expression l|33|) we can 
asymptotically neglect the component in the subspace H- and drop the modulus from \ Y\. In this way, the probability 
distribution ll-iol) can be approximated as 

P^Qdxdre-^i'^^'y 12 dsdre-. (37) 
7r (ia\ Y \ia) 



FIG. 1: Optimal probability distribution for the joint estimation of the squeezing parameter and real displacement for the 
vacuum state. 




FIG. 2: Optimal probability distribution for the joint estimation of the squeezing parameter and real displacement for a coherent 
state with amplitude a = lOi. 



Neglecting the higher order terms, we thus obtain the Gaussian distribution 

p a (x,r)dxdre- r = - e'^^ e~ x * dxdr . (38) 

Notice that, asymptotically the most likely values of the unknown parameters x and r are the true ones x = r = 0, 
and, in addiction, also the mean values of x and r coincide with the true one, namely the estimation is unbiased. 

Prom the asymptotic expression (|38[) can see that the uncertainty in the estimation of the squeezing parameter 
r goes to zero with the number of photons n = a 2 , namely the r.m.s error is Ar = l/\/2fT, while the uncertainty 
in the estimation of the displacement x remains fixed, with the value Aa; = l/\/2- It is interesting to compare the 
precision achieved by the joint estimation with the precision that could be achieved if the parameters x and r were 
measured separately. First, it is known that the optimal estimation of real displacement is given by the observable 
X = (a + a^)/2, and the corresponding uncertainty in a coherent state is given by Ax opt = 1/2. The optimal 
estimation of real squeezing has been recently derived in [2^, and we report here the asymptotic distribution for an 
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excited coherent state: 



2 a 



2 



■2 2 



p(r) = d rm. e' 2 ^ r ~ . (39) 



Accordingly, the r.m.s. error for the optimal estimation of squeezing in a coherent state is Ar opt = l/(2y/n). It is 
remarkable to note the relation 

Ax = V2 Ax°* 

Ar = V2Ar°P* [ ' 

between the uncertainties in the joint measurement and the ones in the separate measurements of squeezing and 
displacement. In particular, the product of the uncertainties in the joint estimation is twice the product of uncertainties 
in the optimal separate measurements: 

AxAr = 2Ax opt Ar opt . (41) 

Surprisingly, this is exactly the same relation as the one occurring in the optimal joint measurement of two conjugated 
quadratures X and Y that can be achieved by heterodyne detection [30j | . 



D. Joint estimation for a displaced squeezed state 



In the previous section we analyzed the optimal estimation of squeezing and displacement for an excited coherent 
state. In that case, while the error in the estimation of r goes to zero with the number of photons, the error in 
estimating x remains fixed. However, it is possible to choose the input state in such a way that both variances vanish 
in the asymptotic limit. To this purpose, we consider here displaced squeezed states \ia, z) — D(ia)S(z)\0) with 
a, z € R. Such states have the wavefunction 



2e 2 



2z\ 1/4 



(y\ia,z)=^—j e"^ e , (42) 

namely a Gaussian centered around the mean value a, with standard deviation a = l/(v / 2e 2 ). Clearly, if the 
conditions «> 1 and a 3> cr are simultaneously satisfied, such a Gaussian lies almost completely in the positive half- 
line. Therefore, in the asymptotic limit a — > +oo, a e~ z , the optimal probability distribution can be approximated 
as 

/ n _ 1 Kjgrg! Ux, r Y \ia,z)\ 2 
7r {ia, z\ Y \ia, z) 

as in the case of coherent states. By calculating the expectation values and keeping the leading order terms, we then 
obtain the asymptotic distribution 

Paz (x,r)dxdre- r = -e^*^ e'^*? dxdr . (44) 

Again, in the asymptotic limit the most likely values in the probability distribution coincide with the true ones, and, 
moreover, the estimation is unbiased. 

The r.m.s. error in the estimation of squeezing and displacement are now given by Ar = l/(V2ae z ) and Ax = 
l/(V2e~ z ), respectively. In order to have both errors vanishing, one needs simultaneously ae z 1 and e~ z ^s> 1. 
For example, we can have an isotropic distribution Ar = Ax with the choice a = e~ 2z . In the isotropic case, we 
notice that only a small fraction of order \[fi of the total number of photons n = a 2 + sinh 2 z comes from squeezing. 
Since one has Ax = Ar = l/(v / 2e _z ) w l/(\/2n 1 / i ), the convergence to the asymptotic regime is quite slow: the 
uncertainty goes to zero with order 1/n 1 ^ 4 . As an example in the asymptotic regime, in Fig. 3 we plot the exact 
probability distribution for n = 4000. 

Also for displaced squeezed states one can compare the accuracy of the joint estimation with that of independent 
measurements of squeezing and displacement. For squeezed states of the form (|42|l the uncertainty in the measurement 
of the observable X is Ax opt = l/(2e~ 2 ). On the other hand, the optimal estimation of squeezing is obtained 
asymptotically by the observable ln(|y|/a), and its uncertainty is Ar opt = l/(2ae z ) (according to the results of [29j|). 
Again, a factor \[2 relates the standard deviations of the marginals in the joint estimation with those of the optimal 
separate measurements. Again, we find the relation (|41|l between the product of uncertainties. 
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FIG. 3: Optimal probability distribution for the joint estimation on a displaced squeezed states with total number of photons 
n = 4000 (with signal photons n — y/E, and squeezing photons \fn) . 



We then see that the estimation of real squeezing and displacement on input states of the form (|42|l corresponds in the 
asymptotic regime to the joint measurement of the observables X and (In | Y \/a), whose commutator is [X, ln(| Y \/a)] — 
i/(2Y) . We thus have the Heisenberg- Robertson inequality 



AXAln(|Y|/a) > 



1 

4Y 



(45) 



where ( ) denotes the expectation value. From this point of view, the displaced squeezed states (|42|l are characterized 
asymptotically as minimum uncertainty states, since they saturate the inequality, and the product of uncertainties 
in the maximum likelihood estimation is exactly twice the Heisenberg limit. To the best of our knowledge, the 
present case is the first example of joint measurement of two harmonic-oscillator noncommuting observables whose 
commutator is not a c- number. 



E. Joint estimation of squeezing, displacement, and reflection 



The results presented in the previous sections can be easily extended to include the estimation of the reflection on 
the real line that is realized by the parity operator P = (— l)" a . In this case we are interested in estimating the three 
parameters in the transformation 

|^) _> P *D(x)S(rM , (46) 

where e can assume the values or 1. The representation {P e D(x)S(r)} is now irreducible in TL, and the associated 
DMC operator is D = 7r/|Y|. For a given input state \ip), the optimal POVM can be written as M(e, x, r) = 
P e D{x)S(r)\r ] }(r ] \S{ry D(xY P e , where 

\V) = , W) , (47) 
yfr ¥\ \Y\ W) 

according to Eq. (|21|) . Clearly, also in this case it is possible to enhance the sensitivity of detection by increasing the 
average value of the modulus |Y|, using coherent states or displaced squeezed states. 



F. Perfect detection of squeezing and displacement: a two-mode setup 

A fundamental mechanism leading to the optimal estimation of group transformations is the use of equivalent 
representations of the group, via the technique of entanglement between representation spaces and multiplicity |l5lll6j . 



12 



Such a strategy strongly improves the quality of estimation for unimodular grou ps JI2T IT3 . ll5T ] , and in the case of 
compact groups is a way to obtain the best estimation of a unitary transformation[l(j. Equivalent representations of 
the affine group can be obtained by entangling the radiation with a second reference mode, which is not affected by 
the unknown squeezing and displacement. This corresponds to considering the two- mode Hilbert space 7i®TL and 
the representation {U XjT <g> 1} of the affine group, namely the reference mode plays the role of an infinite dimensional 
multiplicity space. In this case, the optimal estimation can be obtained with the POVM specified by the general 
formula i|17fl . A remarkable feature of the two- mode setup is that it is possible to have an orthogonal POVM for 
the estimation, namely there exists an ordinary observable on the extended Hilbert space, associated with the joint 
measurement of squeezing and displacement. In fact, by defining the bipartite vectors 

(48) 



(49) 



l*±» = 4= / 00 dy^M \y)\±y) , 
J -00 

we can construct the orthogonal POVM 

M{g) = U x , r ®t ^|*.»«*.|J Ul r ®t. 

The normalization follows straightforwardly from Eq. 0. Moreover, we have the orthogonality relation 

1 f + °° 

<(<fr s | U x . r ® 1 |$V» = - / dy\y\ (y\U Xtr \sy/s') 
1 



7T 



dy \y\ e -r/2 e -^ye- r S (y - y e - r ) 6 S , S , = 6(x)6(r)8 S;S , , (50) 



where we used Eq. (I52|l of the Appendix and the identity S(y — ye r ) = j^r#(r). The vectors in Eq. (|48() are not 
normalizable. However, similarly to the case of the heterodyne operator, where physical states can arbitrarily well 
approximate the unnormalizable eigenstates [2S|. here one can consider, e.g., states of the form 

|*(A)±)) = N X A atQ+btb |*±)) , (51) 

where < A < 1, N\ is a normalization constant, and b^b is the photon number operator for the auxiliary mode b. 
The states |$(A)±)) approaches the optimal vectors H48|) as A — ► 1, with a correspondent increasing to infinity of the 
average energy of the radiation field. 



IV. CONCLUSIONS 



In this paper we presented the joint estimation of real squeezing and displacement of the radiation field from the 
general point of view of group parameter estimation. The combination of squeezing and displacement provides a 
representation of the affine group u ax + 6" , which is the paradigmatic example of a nonunimodular group. To deal 
with the concrete example of squeezing and displacement, we derived in the maximum likelihood approach the optimal 
estimation of a group transformation in the case of nonunimodular groups, providing explicitly the optimal POVM for 
a given input state. In this analysis, some remarkable features of estimation showed up. Firstly, while for unimodular 
groups the maximum likelihood measurements coincide with the usual square-root measurements, for nonunimodular 
groups the SRM are no longer optimal, namely they do not maximize the probability density of detecting the correct 
value. Moreover, for nonunimodular groups one can optimize the estimation strategy in the maximum likelihood 
approach, but even for the optimal POVM the true value is not the one which is most likely to be detected. To 
reduce the discrepancy between the true value and the most likely one a suitable choice of the input states is required. 
Both these features are in general unavoidable, and their origin tracks back to the presence of a positive unbounded 
operator — the Duflo-Moore-Carey operator — in the orthogonality relations for nonunimodular groups. In the problem 
of joint estimating real squeezing and displacement, all the above effects occur. In particular, for coherent input states 
and displaced squeezed states we observed how an increase in the expectation value of \Y\ gives rise to an improvement 
in the quality of estimation, along with a reduction of the discrepancy between the true value and the most likely one. 
In the mentioned cases the probability distributions for joint estimation become asymptotically Gaussian, and the 
r.m.s. errors Aa; and Ar can be easily calculated. Remarkably, the product of uncertainties in the joint estimation is 
exactly twice the product of uncertainties for the optimal separate measurements of squeezing and displacement, in 
the same way as in the joint measurement of two conjugated quadratures. Finally, the use of entanglement with an 
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additional mode of the radiation field allows one to perform a von Neumann measurement for the joint estimation of 
squeezing and displacement, in terms of an ordinary observable with continuous spectrum. The ideal input states for 
detecting an unknown affine transformation can then be approximated by normalizable states. 
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V. APPENDIX 

Here we derive the group theoretical structure of the representation {U x . r = D(x)S(r)} of the affine group, by 
explicitly calculating the expression for the group average of an operator A over the left-invariant measure d^ g = 
dxdre 



— r 



Proposition 3. The Clebsch- Gordan series of {U x , r } contains two irreducible irreps, + and — . Accordingly, the 
Hilbert space can be decomposed as TL — TL+ © H.-, and the projections onto the irreducible subspaces are 1+ = 0(Y) 
and 1- = 6(—Y), respectively. The DMC operators are given by D+ = tt9(Y) /\Y\ and D- = tt8(— Y)/\Y\. 

Proof. Using twice the resolution of the identity in terms of the eigenstates \y) of the quadrature operator Y, and 
the relation 

U x<r \y) = D(x)S(r)\y) = D{x)e^' 2 \e- r y) = e^e' 2 ^ \e~ r y) , (52) 
we can calculate the group average as follows 

poo pOO 

(A) = I dre~ r / dxU x , r AUl r 

> J — OO 

pOQ pOO pOO 

dre' r / dx / dy / d y' (y\ A\y') e -^iy-v')^ r e - r \ e - r y) (e~ r y'\ 



— oo *i —oo 

OO coo 



tt dr dy(y\A\y)e r \e r y)(e r y\ 



sgn(y)oo 



tt I dy I dr- (y\A\y) \f)(f\ 

/o y 



— oo 



poo />0 

= tt / df \f)(f\ Tr[A0(Y)/Y] - tt / df\f)(f\ Tr[A6(-Y)/Y] 

JO J -oo 

= 6(Y) Tr[A tt9(Y)/\Y\] + 6(-Y) Ti[A tt9{-Y)/\Y\] . (53) 
The thesis follows by comparing the last equation with the general formula (jSJ) for the group average. 
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